Numerical investigation of heat and mass transfer in three-dimensional MHD nanoliquid flow with inclined magnetization

Heat and mass transfer rate by using nanofluids is a fundamental aspect of numerous industrial processes. Its importance extends to energy efficiency, product quality, safety, and environmental responsibility, making it a key consideration for industries seeking to improve their operations, reduce costs, and meet regulatory requirements. So, the principal objective of this research is to analyze the heat and mass transfer rate for three-dimensional magneto hydrodynamic nanoliquid movement with thermal radiation and chemical reaction over the dual stretchable surface in the existence of an inclined magnetization, and viscous dissipation. The flow is rotating with constant angular speed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\upomega }^{*}$$\end{document}ω∗ about the axis of rotation because such flows occur in the chemical processing industry and the governing equations of motion, energy, and concentration are changed to ODEs by transformation. The complex and highly nonlinear nature of these equations makes them impractical to solve analytically so tackled numerically at MATLAB. The obtained numerical results are validated with literature and presented through graphs and tables. Increasing the Eckert number from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$5\le Ec\le 10,$$\end{document}5≤Ec≤10, a higher Nusselt and Sherwood number was noted for the hybrid nanofluid. By changing the angle of inclination \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α, the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Nu}_{x}$$\end{document}Nux performance is noted at 8% for nanofluid and 33% for hybrid nanofluid. At the same time, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${Sh}_{x}$$\end{document}Shx performance of 0.5% and 2.0% are observed respectively. Additionally, as the angle of inclination increases the skin friction decreases and the chemical reaction rate increases the mass transmission rate.

In a myriad of industrial applications, ranging from nuclear reactors and automobiles to electronics, the role of fluids is pivotal in enhancing heat transfer rates (HTR).Traditional fluid options like water, oils, and ethylene glycol, however, exhibit limited thermal conductivity.This limitation has spurred extensive research endeavors aimed at elevating HTR for improved efficiency and performance.First, Choi and Eastman 1 coined the term "nanofluid (NF)" by mixing non-metallic or metallic nanoparticles (NPs) in host fluid which dramatically augments HTR.The characteristics of such fluids depend upon different factors like the shape and size of suspended nanoparticles in host fluid, thermal conduction, etc.Later, many investigators [2][3][4] followed his idea and worked on the enhancement of HTR.Similarly, the mixing of two or more nano-sized particles in the host fluid is known as a hybrid nanofluid (HNF).Continuous and different strategies have been adopted to enhance HTR.Jena et al. 5 explored the recent development in heat transfer (HT) characteristics of NF by considering the non-uniform heat source and inclined magnetization.Parida et al. 6 computationally discussed the dust particles in water and kerosene-based nanoliquid for HTR.Pattnaik et al. 7,8 used different NPs like copper, aluminum, gold, and singlewall carbon nanotubes to explain HTR by using water as host fluid over permeable surfaces.
In industrial sectors like aerodynamics, plastic sheet extrusion, continuous metallic plate extrusion, artificial fiber synthesis, and plastic film magnification, the motion of an incompressible fluid over an expanding surface is a frequently observed occurrence.The HTR at the deformable surface plays a substantial role in determining the overall quality of the product in each of these applications.Sakiadis 9 gave the thought of boundary layer (BL) enhancement during the movement of the surface and attracted researchers' attention.The Sakiadis' problem was expanded by Erickson et al. 10 who also looked at the causes of puffing or sucking at the moving sheet on HMT in the BL flow.Baag et al. 11 for exploration of MHD boundary layer flow over porous exponentially SS with a uniform heat source.Nayak et al. 12 considered the radially stretched sheet for discussion by incorporating the variable magnetic field.Mishra et al. 13 investigated the HT influence on the MHD movement of micropolar liquid passing from a permeable medium considering a similar heat source.Upreti et al. 14 used the Casson NF over SS utilizing the Cattaneo-Christov model in stagnation point flow to examine the shape factor.Seini and Makinde 15 explored the magnetohydrodynamic (MHD) BL flow above the exponentially stretched sheet by considering the chemical reaction (CR).Arshad et al. 16 considered radiative heat and mass transfer (HMT) for HNF incorporating the inclined magnetic field.Upreti et al. 17 explored the influence of shape factor on Casson goldblood nanofluid flow through SS incorporating magnetic effect.Also, considered the impact of Ohmic heating, convective heating, and suction/injection on heat transfer rate.Singh et al. 18 investigated the influence of Melting and CR on the immobility point flow of a micropolar fluid over a Porous SS Medium.Pandey et al. 19,20 presented the multiple slip mechanism and volumetric heat generation over porous SS and cone respectively.Sreedevi and Reddy 21 considered 3D NF flow above SS with radiation and CR for the exploration of thermo-diffusion and Brownian movement.They concluded that enhancing the Deborah number increases the temperature profile.
In recent years, significant research attention has been directed toward Heat and Mass Transfer flows due to their essential nature, prevalent in various engineering and industrial sectors.An example is the extrudate from the die in a melt spinning process.Researchers have investigated multiple flow scenarios, including those involving a stretching surface (SS), chemical reaction (CR), mixed convection, and thermal radiation, to address HMT-related issues.Rao et al. 22 delved into HMT aspects in the context of a thermally oscillating fluid.Arshad et al. 23,24 explored HMT by considering thermal radiation and chemical reactions for different types of Newtonian fluids (NFs) and Hybrid Nanofluids (HNFs).Mathur et al. 25 examined the Darcy-Forchheimer skin coefficients and explored the velocity slip properties of a micropolar NF.Upreti et al. 26 explored thermodynamics and HT using the Riga plate for the magnetized Casson HNF.They also discussed the entropy generation.Jayavel et al. 27 provided a discussion on heat transfer analysis and irreversibility in MHD Darcy-Forchheimer movement of Casson HNF flow over wedge and cone.They used nanofluid and hybrid nanofluid for their discussion.Hassan et al. 28 used the molybdenum di-sulfide NPs to explore the HMT incorporating the non-linear and linear radiation.Hussain et al. 29 computationally investigated the thermal radiation to find the HTR over a stretchy surface.Different related research [30][31][32][33][34][35][36][37][38] in literature can be found.Arshad and Hassan 39 studied the heat and mass transmission rate in a rotating permeable system using hybrid nanofluids.
Nowadays, chemical reactions, thermal radiation, and the presence of heat source/sink are fundamental components in the study of HMT phenomena.In various engineering applications, the interplay of chemical reactions, thermal radiation, and heat source/sink mechanisms plays a pivotal role in shaping the thermal behavior of systems.Understanding the effects of chemical reactions and thermal radiation, along with heat source/ sink interactions, is essential for optimizing processes in fields such as materials science and chemical engineering.Heidary et al. 40 numerically investigated the magnetic field effect with forced convection in a duct for NF flow.Sheikholeslami and Rokni 41 gave a review on the simulation for the HT phenomenon of nanoliquid in the existence of a magnetic field.Makinde and Mishra 42 examined the MHD mixed convection with non-uniform viscosity Blasius flow inserted in a permeable medium incorporating the chemical reaction.Using the CR and heat source, a semi-analytical solution for MHD Jeffery fluid flow is provided by Nisar et al. 43 .Mishra et al. 44 explained the influence of nonlinear radiation and cross-diffusion effects on the flow of micropolar nanoliquid over a stretching sheet with an exponential heat source.Mehrizi et al. 45 reported a new analysis of natural convection BL flow with variable wall temperature on a horizontal plate.Reddy et al. 46 explored the HMT flow of NF at inclined plates with thermal radiation and magnetic fields under enhanced boundary conditions.Nayak et al. 47 investigated the flow and HTR non-Newtonian fluid with hybrid nanoparticles by employing a magnetic field.
The literature review conducted indicates a notable gap in research, as there has been no investigation into comparing different types of nanofluids over dual stretchable surfaces while factoring in the existence of an inclined magnetic field and accounting for viscous dissipation.Such applications involving nanofluids occur in advanced cooling systems, enhanced heat exchangers, biomedical devices, material processing, renewable energy, etc. Jena et al. 5 considered only the temperature profile in their study but ignored the concentration profile.Jayavel et al. 27 made their analysis over the wedge and cone but ignored the inclined magnetization.Pattanaik et al. 7 examined the uniform heat source and ignored the chemical reaction.Gupta et al. 32 used the kerosene oil for their investigations over exponentially SS.Arshad and Hassan 39 deliberated hybrid nanofluids using different NPs.Singh et al. 18 considered the non-uniform heat source but ignored the thermal radiation.Similarly, Pattanaik et al. 38 , Mohanty et al. 31 , and Parida 6 presented their research in different aspects but none of these considered the inclined magnetic field.Based on the literature survey conducted above, the novelty of this study aims to provide a comparative analysis of different water-based NFs.These NFs, namely Cu/H 2 O ; NF, andCu − Al 2 O 3 /H 2 O HNF, are investigated concerning their HMT rates.This report is prepared for a dual stretching surface placed within a porous medium.Viscous dissipation, thermal radiation, chemical reaction, and most importantly inclined magnetic field are considered because they involve semiconductor manufacturing, solar energy, food processing, etc.The governing equalities are changed into the ODEs by employing a transformation and tackled a MATLAB by BVP-4c algorithm by setting the tolerance 10 −6 for solutions.The results are obtained for increasing values of different parameters involved in this research.This comparative analysis serves to provide insights into the central research queries outlined below: 1. How does the rotation parameter affect the velocity, temperature, and concentration profile? 2. What is the influence of increasing behavior of mixed convection, stretching ratio, magnetic force, and inclination angle on velocity and temperature profile? 3. Does the increased thermal radiation and viscous dissipation reduce skin friction and enhance the HTR? 4. How does the chemical reaction affect the concentration profile and Sherwood number? 5. What do we get numerical outcomes for skin frictions along thex − axis , y − axis, Nusselt, and Sherwood numbers versus different parameters?

Problem formulation
Consider a steady three-dimensional boundary layer, the MHD flow of a fluid over a dual stretching sheet within a porous medium, accompanied by a chemical reaction, thermal radiation, and the existence of a nonuniform heat source.The conceptual representation of the problem, inclusive of the flow design and coordinate system, is presented in Fig. 1.The horizontal direction is represented by the x-axis, the upward direction by the z-axis, and the y-axis is perpendicular to both axes.The fluid is undergoing a steady rotation at a constant speed denoted as ω * around the z-axis .The surface is stretching along the x-axis with velocity U w = ax and along the y-axis with velocity V w = by as shown.Inclined magnetic field B 0 with angle α from the x-axis to the axis of rotation is working.By these assumptions, the flow conservation, momentum along the x-axis, and y-axis, temperature and concentration in the existence of joule heating, thermal radiation, and viscous dissipation equation are: Equation of continuity 48 : www.nature.com/scientificreports/Momentum equations along the x and y axis 48 : Energy equation without q r relation 48,49 : Concentration equation 49 : The respective boundary conditions 48 for the current problem are: Here g * -gravitational acceleration, B 0 -magnetic field, T-temperature, C-concentration, T ∞ -ambient temperature, q r -radiative heat flux, µ-dynamic viscosity, ρ-density, D B -mass diffusion, D T -temperature diffusion, k-thermal conductivity, C p -specific heat, σ-electrical conductivity, B t -thermal volumetric coef- ficient, α-the angle of inclination, u, v, w are velocity components in x, y, z respectively, the subscript hnf rep- resents the hybrid nanofluid.
The specified boundary conditions for the fluid flow and thermal transport problem outline the physical behavior near and far from the solid surface.At the origin (z = 0), the prescribed velocities U w and V w denote a stretching wall, while w = 0 enforces a no stretching.Temperature (T) and concentration (C) conditions T w and C w at the wall capture heat and mass transfer interactions.As z → ∞, the velocity components approach zero (u → 0, v → 0), signifying a quiescent state, while temperature and concentration www.nature.com/scientificreports/ the free-stream values far from the surface.These conditions collectively provide interactions and asymptotic behavior of fluid flow, heat transfer, and mass transport in the given problem.
There are numerous models in literature for characterizing the effective properties of NF and HNF.So, the thermophysical relations of NF and HNF are given in Tables 1 and 2 correspondingly.Table 3 shows the thermophysical values for used NPs and base fluid.Table 4 shows the comparison of literature and present outcomes.

Properties
Hybrid nanofluid relations 39 Density Here Electrical conductivity   The term q r in the right side of Eq. ( 4) presents the heat radiation effect.Utilizing the Rosseland approxima- tion, the radiative heat flux is computed in the following: Here σ * represents the Stefan-Boltzmann coefficient and k 1 represents the mean absorption constant.The radiation is optically thick considered.Considering that sufficiently small temperature difference in the flow, the term T 4 by employing the Taylor series as follows: Therefore, by abandoning higher-order terms above the first degree in (T − T ∞ ), we get.By using Eqs.( 7) and ( 8) So, Eq. ( 4) takes the following form: Similarity transformation: In the present problem, the following similarity transformations are taken to change the PDEs into ODEs.
The differentiation is w.r.t η.The flow Eq. (1) of mass conversation is satisfied identically by employing Eq. ( 12).The Eqs. (2,3,5) and Eq. ( 11) for momentum, energy, and concentration will take the following form after using Eq. ( 12) with boundary conditions: The non-dimensional quantities H 1 , H 2 , H 3 , H 4 , H 5 , and K 2 , are hybrid NPs relations (presented in Table 2), and , δ, Z, ǫ, M, Pr, π, Ec, Sc, N t , τ , N b , and K c are defined as (T − T ∞ )x 3 , Re x = u w v f .The modified boundary conditions are presented as follows: Here γ = b a presents the dimensionless stretching ratio.

Engineering parameters of interest:
There are the following most important quantities regarding to engineering perspective.(A) Skin friction: The significant surface skin friction along the x-axis and y-axis are Cf x , Cf y defined as: The τ zx and τ zy indicate shear stress along the stretched wall x-axis and y-axis which is defined as: The dimensionless form of Eq. ( 20) with the help of Eq. ( 16) is: (B) Rate of heat and mass transfer: The HMT rates, expressed through Nusselt and Sherwood numbers are the following: By utilizing the temperature field, the thermal diffusion rate is characterized by the Nusselt number: By using the similarity transformation, the transformed form of Eq. ( 22) is given below: By utilizing the concentration field, the mass transmission rate is characterized by the Sherwood number: By using the similarity transformation, the transformed form of Eq. ( 24) is given below:

Numerical solution
This section provides the numerical solution methodology as presented in Fig. 2 to get the solution of higherorder nondimensional ODEs.For this purpose, the bvp4c method is used to tackle the boundary value problem by transforming it into an initial value problem.The Fig. 3 validates the code of the present problem.For this purpose, newly defined variables are as follows: The following form of equations is used in MATLAB to get the numerical solution.

Results and discussion
The problem is formulated for HNF in the previous section.The authors obtained the results of the current problem separately for both NF and HNF and discussed them physically.The obtained influences of different constitutes by using the bvp4c algorithm at MATLAB are described in the following.

Velocity profiles
The influences of study parameters on velocities are discussed in this section.The influence of rotation on velocity profile p′(η) when is presented in Fig. 4a.The graph displays the velocity profiles p′(η), demonstrating a decrease in veloc- ity as the rotation parameter increases for both NF and HNF cases.Initially, at = 0 and = 1, there is minimal alteration in the velocity profiles.However, a significant acceleration in the decay of these profiles is observed as the rotation parameter rises to = 2 and = 3.The main reason behind this is enhanced viscous dissipation that affects the velocity profile.Figure 4b represents the influence of the magnetic field on velocity q′(η)whenNb = 0.1, The magnetic parameter has a great influence on the velocity profile.In the non-existence (M = 0) of the magnetic field, the momentum BL looks similar but when it upsurges between 1 ≤ M ≤ 3, a big change is noted due to the existing Lorentz force which decreases the fluid motion and consequently momentum BL decay.A wider momentum BL is noted for HNF as associated with NF.This is due to the generation of damping electromagnetic force which restricts the fluid's motion and results in a gradual reduction in velocity profile.Figure 4c   , M = 3 and a higher velocity profile is observed for NF as associated with the HNF.This is owing to the driving forces associated with convection which typically increases the fluid's velocity.Figure 4e shows that porous medium parameter Z has a negligible decreasing influence on velocity profile p′(η) . The porosity acts as a flow impediment, and the drag force significantly reduces the velocity of the fluid.The stretching ratio parameter γ represents the effects on velocity profile q′(η) in Fig. 4f  whenNb When the stretch- ing rate along the x − axis increases as compared to the y-axis , the velocity momentum BL expands for both types of NFs.

Temperature and concentration profiles
The following Fig. 5a-h show the impression of different constraints on temperature and concentration profiles.The influence of magnetic parameters on temperature r(η) is represented in Fig. 5a  whenNb seen that the temperature increases when magnetic strength increases.The justification for this is that Lorentz force reduces the fluid motion and as a result, a maximum HT and a higher thermal BL are noted for HNF as compared to NF. Figure 5b indicates the effect of rotation on the temperature (η) whenNb = 0.1, δ = 0.8, M = 0.3, Z = ǫ = π = Kc = Ec = 0.5, Nt = Le = = 0.3, α = 45 • .W h e n t h e rotation of fluids increases the temperature profile shows the same behavior.Rotation increases temperature transmission by promoting fluid mixing, reducing boundary layer thickness, and enhancing convective heat transfer due to the formation of vortices and turbulence, leading to improved thermal.In this figure low but consistent HTR is noted for HNF.The influence of thermal radiation on the temperature r(η) is shown in Fig. 5c whenNb In the absence of the radiation parameter, a higher temperature BL is detected for both NFs.Thermal (34) radiation typically decays temperature profiles by emitting heat energy from hotter regions to cooler surroundings, causing temperature gradients to diminish, resulting in a more uniform temperature distribution.The Eckert number characterizes the thermal energy conversion associated with the fluid's motion.A higher Eckert number implies more kinetic energy, which promotes cooling, leading to a decrease in temperature gradients and a more uniform temperature profile.This phenomenon is presented in

Numerical outcomes
The numerical outputs for skin frictions Cf x , Cf y , Nusselt number Nu x , and Sherwood number Sh x are presented in this section.Table 5 represents the obtained numerical results for skin frictions Cf x and Cf y for both nanoflu- ids.A minimum skin friction Cf x is observed in the absence of rotation i.e. λ = 0 and increasing behavior is seen for increasing values of rotation parameter from = 1 to = 2, while an opposite behavior is observed for Cf y .
Increased porosity from 10 ≤ Z ≤ 30 reduces skin friction in porous media because higher permeability allows fluid to flow more easily through the medium, reducing the resistance to motion and minimizing the friction between the fluid and the porous structure.The stretching ratio is defined by the rate along the y-axis to the  www.nature.com/scientificreports/stretching rate x-axis .Increased stretching ratio 0.9 ≤ γ ≤ 1.1 increases skin friction in the context of flow over a solid surface.This occurs because a higher stretching ratio leads to increased shear stress.Reduced skin friction is followed in the deprivation of mixed convection and magnetic field and increasing behavior is noted when it increases.The angle of inclination has a great influence on skin friction.When the magnetic field is parallel to the x-axis , a maximum skin friction Cf x is noted while it decays when the angle of inclination increases between the range 0 ≤ α ≤ 90.Minimum skin friction is noted when the magnetic field becomes parallel to the y − axis.
Table 6 shows the influences of different constraints on Nusselt number Nu x and Sherwood number Sh x for both nanofluids.A higher HMT rate is noted when the angle of inclination for the magnetic field is parallel to the xaxis.Increasing radiation parameter 0 ≤ π ≤ 4 typically decreases the Nusselt and Sherwood numbers because radiative heat or mass transfer is a slower process compared to convection.As radiation becomes more significant, it diminishes the relative contribution of convection, leading to lower Nusselt and Sherwood numbers, indicating reduced HMTR and utmost Nusselt and Sherwood is seen in the absence of radiation, i.e. π = 0 .Prandtl number 7.2 ≤ Pr ≤ 9.

Conclusion
In this present article, the authors discussed the heat and mass transfer rate of 3D rotating, incompressible NF and HNF fluid over the dual stretchable surface.The governing equations are transformed and tackled in MAT-LAB.We considered the inclined magnetic field, viscous dissipation, thermal radiation, and chemical reaction

Figure 1 .
Figure 1.Flow configuration of the problem.

Figure 2 .
Figure 2. Flow chart for the numerical solution.

2
and field 0 ≤ M ≤ 15 have a direct relation and inverse relation with HMTR respectively because Pr implies that thermal or mass diffusion is rapid relative to momentum diffusion and the existence of magnetic field suppresses fluid movement, disrupting convection resulting in reducing HMTR.Increasing values of the CR parameter 0 ≤ Kc ≤ 4 and Lewis number 10 ≤ Le ≤ 20 decreases the Nusselt number while it increases the Sherwood number.CR consumes the temperature and concentration gradient and a higher Lewis number results in a more dominant effect on fluid behavior.Similar results are found for thermophoresis 0 ≤ N t ≤ 6 and Brownian motion 0.5 ≤ N b ≤ 1.5 factors corresponding to the Nusselt number and Sherwood number.As the viscous dissipation parameter increases between 5 ≤ Ec ≤ 15 the HMT rate increases in both cases of used nanofluids because increased Ec number indicates that the kinetic energy of the fluid is relatively higher than internal energy change.

Table 4 .
Comparison of present outcomes with literature.